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ABSTRACT 

Aims. We present a 8° x 6°, high resolution extinction map of the Pipe nebula using 4.5 million stars from the Two Micron All 
Sky Survey (2MASS) point source catalog. 

Methods. The use of Nicer (Lombardi & Alves 2001), a robust and optimal technique to map the dust column density, allows 
us to detect a Av = 0.5 mag extinction at a 3-a level with a 1 arcmin resolution. 

Results, (i) We find for the Pipe nebula a normal reddening law, E{J — H) = {l.S5±0.15)E{H ~ K). (li) We measure the cloud 
distance using Hipparchos and Tycho parallaxes, and obtain I30I5I P*^- This, together with the total estimated mass, 10"* M©, 
makes the Pipe the closest massive cloud complex to Earth, (lii) We compare the Nicer extinction map to the NANTEN 
^■^CO observations and derive with unprecedented accuracy the relationship between the near-infrared extinction and the ^'^CO 
column density and hence (indirectly) the ^^CO X-factor, that we estimate to be 2.91 x 10^" cm~^ km~^ s in the range 
Av £ [0.9, 5.4] mag. (iv) We identify approximately 1 500 OH/IR stars located within the Galactic bulge in the direction of the 
Pipe field. This represents a significant increase of the known numbers of such stars in the Galaxy. 

Conclusions. Our analysis confirms the power and simplicity of the color excess technique to study molecular clouds. The 
comparison with the NANTEN ^'^CO data corroborates the insensitivity of CO observations to low column densities (up to 
approximately 2 mag in Av), and shows also an irreducible uncertainty in the dust-CO correlation of about 1 mag of visual 
extinction. 

Key words. ISM: clouds, dust, extinction, ISM: structure, ISM: individual objects: Pipe molecular complex. Methods: data 
analysis 



1. Introduction 

Nearby Galactic molecular clouds complexes represent our 
best chance to understand cloud formation and evolution 
and hence to learn how stars come to be. But progress 
on the study of these objects has been slow. Not only 
are molecular clouds the coldest objects known in the 
Universe, their main mass component (H2) cannot be de- 
tected directly. Most of all we know today about their 
physical properties has been derived through radio spec- 
troscopy of H2 surroga t e molecules (CO, CS, NH3; e.g., 
iBlitz fc Winiamd Il999l: iMyerd Il999(l and more recently 
through thermal emission of t he dust g rains inside these 
clouds ( Andr e et a l. 2000; Jo hnstone et a l. 2000). The re- 
sults obtained using these techniques, especially the esti- 
mate of column densities, are not always straightforward 
to interpret and are plagued by several poorly constrained 
effects. Moreover, although large scale maps of entire 
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mole c ular cloud complex es are now available IjHever et alJ 
Il998t ISimon et al.ll200lh . maps with sufficient resolution 
and dynamic range to identify not only dense molecular 
cores but also their extended environment are still not ex- 
istent. This wide view on molecular clouds is at present an 
obvious gap in our understanding of the relation between 
the dense Interstellar Medium (ISM) and star formation. 

A straightforward and powerful techniq ue to study 
molecular cloud structures, pioneered by iLada et "all 
l|l994j) and kn own as the Near -Infrared Color Excess 
method (Nice, lAlves et al.lll99?i) makes use of the most 
reliable tracer of H2 in these clouds: extinction by perva- 
sive dust grains in the gas. The Nice method relies on 
near-infrared (NIR) measurements of extinguished back- 
ground starlight to derive accurate line-of-sight estimates 
of dust column densities. Depending on the stellar rich- 
ness and color properties of the background field this tech- 
nique can produce column densi ty maps with spat i al res- 
olutions down to 5 arcsec (e.g. lAlveset allboOll l2002t) 
and with dynamic ranges more then an order of magni- 
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tude larger than classical optical star count techniques. 
This novel view on molecular clouds is providing not only 
new infor mation on the physical structure of these object 
ijAlves et al. 2001') but also an insight into their chemical 
structure and the physical properties of the dust grains, 



when combined to molecular line and dust emission data 
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In part to make use of the wealth of NIR data 
provided by the Two M icron All Sky Survey (2MASS; 
iKleinmann et alJ Il994|) the Nice method was fur- 
ther developed into an optimized multi-band technique 
dubbed Near-Infrared C olor Excess Revisited (Nicer, 
iLombardi fc AlveslEoOll hereafter Paper I). This gener- 
alization of Nice can, nevertheless, be applied to any 
multi-band survey of molecular clouds. Through use of 
optimal combinations of colors. Nicer improves the noise 
variance of a map by a factor of two when compared to 
Nice. This unique property of Nicer makes it the ideal 
tool to trace large scale distributions of low column den- 
sity molecular cloud material. When applied to 2MASS 
data, Nicer dust extinction maps trace not only the low 
column density regions {Ay ~ 0.5 mag) but have the dy- 
namic range to identify dense molecular cores by reaching 
cloud depths of ^ 30 m ag of extinction, c orresponding to 
8 X 10^2 protons cm-^ |)Lombardill2005hj) . 

In this paper we present an extinction map of the 
Pipe nebula covering 48 sq deg, computed by applying the 
Nicer technique on 4.5 millions JHK photometric mea- 
surements of stars from the 2MASS database. The Pipe 
nebula is a poorly studied nearby complex of molecular 
cloud s. The only systema tic analysis of this region is the 
one of lOnishi et alJ l)l999() who present a ^ 27 sq deg map 
in the J = 1 — line of ^^CO observed on a 4 arcmin grid, 
and smaller maps of selected regions in the J = 1 — lines 
of ^^CO and C^^O. They estimate the total ^^CO mass to 
be ~ 10"* Mq (for a cloud distance of 160 pc) and point 
out that star formation is only occurring on Barnard 59, 
where the C-'^^O column density is the highest and where 
they find a CO ou tflow. Bar n ard 59 was also observed 
at 1300 /xm bv iReipurth et al.l l)l996(l and a protostellar 
candidate B59-MMS1 was found. Two pre-main sequence 
stars associated with Barnard 59 appea r in the young bi- 
naries sur vey of IReipurth fc Zinnecked l|l993j) and more 
recently in iKoreskd l)2002j) " 

This paper is organized as follows. In Sect. |5] we de- 
scribe the technique used to map the dust in the Pipe neb- 
ula and we present the main results obtained. In Sect. |3 
we address the determination of the cloud distance us- 
ing Hipparcos data. A statistical analysis and a discus- 
sion of the bias introduced by foreground stars is pre- 
sented in Sect.lH W e compare the CO observations from 
lOnishi et al.l l)l999(l with our outcome in Sect.|Sl Sectional 
is devoted to the mass estimate of the cloud complex. 
Finally, we summarize the conclusions of this paper in 
Sect.EI 




Fig. 1. Color optical image of the region around the Pipe 
nebula. The bright "star" on the top-right is Jupiter; be- 
low it are well visible the p-Ophiuci dark cloud and the 
Lupus complex. The white box encloses the field stud- 
ied in this paper. The Pipe nebula obscures part of the 
Galactic bulge, and thus occupies an exceptional loca- 
tion for infrared studies (image courtesy Jerry Lodriguss, 
http : //www. astropix. com/). 

In the following we will normally express column densi- 
ties in terms of the 2MASS band extinction Ak ■ When 
converting this quantity into the widely used visua l extinc- 
tion Av, we will use the lRieke fc Lebofskvl l|l985|) redden- 
ing law converted in to the 2MASS photometric system 
^see lCarnent'CTll200ll) . Ak/Av ^ 0.112. 

2. Nicer extinction map 

The data analysis was carried out using the Nicer method 
described in Paper I (to which we refer for more detailed 
information). Infrared J (1.25 /^m), H (1.65 /^m), and 
Kg band (2.17 fim) magnitudes of stars in the Pipe re- 
gion were obtained from the Two Micron All Sky Survey^ 
('2MASS: IKleinmann et"ai]ll994|) . We selected a large re- 
gion around the Pipe nebula, characterized by galactic 
coordinates 

-4° < ; < +4° , +2° < 6 < +8° . (1) 

As shown by these coordinate ranges, the Pipe nebula oc- 
cupies an exceptionally good location for infrared studies. 
Indeed, this cloud complex is just above the galactic plane 
and is in front of the galactic bulge (see Fig.^. The cloud, 
we anticipate, is approximately at 130 pc (see Sect. 131, and 
thus almost all the stars observed in direction of the cloud 
are background objects. As a result, we could carry out 
the study of the infrared extinction of the cloud with neg- 
ligible contamination by foreground objects (except in the 
most dense regions). 

The high density of stars in the Pipe region has also 
some drawbacks. Indeed, because of confusion, the nom- 
inal 2MASS photometric completeness limits drop close 

^ See http://www.ipac.caltech.edu/2mass/. 
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Fig. 2. Color-color diagram of the stars in the Pipe nebula 
field, as a density plot. The contours are logarithmically 
spaced, i.e. each contours represents a density ten times 
larger than the enclosing contour; the outer contour de- 
tects single stars and clearly shows a bifurcation at large 
color-excesses. 

to the galactic center. In the case of Nicer (in contrast 
with the star-counting method), this does not affect the 
measurements of column densities because no assumptions 
are made about the local star density. High density regions 
also pose severe challenges for data reduction. For exam- 
ple, the pipeline used for the Second Incremental Release 
of 2MASS produced slightly inaccurate zero-points. This 
problem showed on the Nicer maps of the Pipe nebula 
as abrupt artificial changes in the measured column den- 
sity at the boundaries of the 2MASS observation tiles. 
Fortunately, the pipeline has been greatly improved and 
this problem does not appear in the final 2MASS All Sky 
Release. 

After selecting point sources from the 2MASS cata- 
log inside the boundaries we generated a preliminary 
extinction map. As described in Paper I, this map was 
mainly used as a first check of the data, to select a control 
region on the field, and to obtain there the photomet- 
ric parameters to be used in the final map (see Fig. ^ . 
We note that the search of a control field close to the 
Pipe nebula has been non-trivial, because of the complex 
cloud structure and high column densities observed at low 
galactic latitudes. We identified in the Eastern part of our 
field (top-right in Fig. O a small region that is appar- 
ently affected by only a negligib le extinction (se e below); 
note that the newly determined llndebetouw et' al. (2005) 
2MASS reddening law was used at this stage. 

Using the information provided by the control field, we 
generated a second map, which is thus "calibrated" (i.e., 
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Fig. 3. Spatial distribution of the subsample of sources 
as defined by Eqs. (O and ©. Subsample A is shown as 
filled circles, while subsample B is shown as open circles 
(see also Fig.|2l. Subsample A appears to be strongly clus- 
tered in high-column density regions of the cloud, and are 
thus interpreted as genuine reddened stars; subsample B 
seems not to be associated with the cloud, and are in- 
stead preferentially located at low galactic latitudes. The 
contour fine represents the Ak = 0.6 contour of the Pipe 
nebula. 

provides already, for each position in our field, a reliable 
estimate of the column density). We then considered the 
color-color diagram for the stars in the catalog. The result, 
shown in Fig.|21 shows a surprising bifurcation for H—K > 
1 mag, which in principle might represents a problem in 
color-excess studies. 

In order to further investigate the origin of this bifur- 
cation, we considered the region H — K > 1 mag in the 
color-color diagram, and divided it into two subsamples 
A (upper branch) and B (lower branch) according to the 
expressions 

A={lA{H~K) + 0.5mag<{J-H)} , (2) 
B = {1A{H -K) + 0.5 mag > ( J - H)} . (3) 

We then analysed the spatial distribution of the stars con- 
tained in the two regions. The results, shown in Fig. 13 give 
a first strong indication that, as expected from their loca- 
tion in the color-color diagram, subsample A is associated 
with the densest regions of the molecular cloud and thus 
it has to be interpreted as normal stars reddened by the 
dust, while subsample B does not appear to be associated 
with the cloud. Interestingly, subsample B is uniformly 
distributed in galactic longitude, but has a strong prefer- 
ence for low galactic latitude regions. 

Further hints on the nature of the two stellar popula- 
tions in subsamples A and B are given by the histogram of 
their K band magnitudes, shown in Fig. 01 Interestingly, 
while subsample A stars show, as expected, a broad distri- 
bution, with number counts increasing at relatively large 
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Fig. 4. The histogram of the K band magnitude for the 
two star subsets A and B of Eqs. (O and 0. 
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Fig. 5. The extinction-corrected color-color diagram for 
the Pipe nebula. The stripes show both subsamples Bi 
and B2 C Bi [of. Eqs. ijU and ©]. 



magnitudes, the B stars show a well defined distribution, 
with a strong peak at ii" ~ 7 mag. 

The lack of correlation between the dust reddening 
and the stars of subsample B is also well expressed by 
extinction-corrected color-color diagram, shown in Fig. El 
This plot has been constructed by estimating, for each 
star of Fig. |21 its intrinsic color, obtained by correcting 
the observed one by the measured Nicer extinction at the 
location of the star (cf. Fig.[7J|. This analysis, based on the 
simplifying assumptions that all stars are background to 
the molecular cloud and that the extinction is basically 
constant on the resolution of our map, is however able to 



capture the essential characteristics of the distribution of 
intrinsic star colors. By comparing Fig. [S] with Fig. [3 in 
particular, we see that the upper branch, subsample A, es- 
sentially disappears in the extinction-corrected plot, while 
the lower branch, subsample B, is largely left unchanged. 
Because of the much shrunken distribution of stars. Fig. El 
is particularly useful to better identify stars belonging to 
the lower branch. We defined thus in this plot two sets 

Si = {\IA{H - X)intr - (J - if)intr | < 0.6 

and {H - K)intr > 0.6} , (4) 

^2= {|l.4(i/-X)intr-(J-if)intr| < 0.6 

and (J - i/)intr > -1-6{H - X)intr + 2.4} . (5) 

Note that B2 C Bi. These two sets, marked in Fig. \S\ 
correspond to the areas in the color-color plot where a 
contamination by subsample B stars is possible (-Bi), or 
highly likely (i?2)- 

From the elements considered so far, we can carry out 
the following conclusions with respect to the nature of the 
B subsample that from this point on will be called the 
"lower branch" : 

— These stars are unrelated to the molecular cloud red- 
dening and appear to be located close to the galac- 
tic plane. A deeper analysis, carried out on a much 
larger field (Lombardi et al., in prep.), actually shows 
that the stars are preferentially located on the galactic 
bulge. 

— The stars occupy a well defined region in the 
(extinction-corrected) color-color diagram. This indi- 
cates that the "lower branch" represents a well-defined 
population of stars. 

— The K band distribution of "lower branch" stars is 
strongly peaked at if ~ 7 mag. This not only is a 
further indication that we are looking at a homoge- 
neous population, but also suggests that these stars 
are at essentially the same distance. The _?C-band stan- 
dard deviation is approximately 0.8 mag (if the tail at 
K > 10 mag of Fig.0]is neglected), and since this must 
account for the scatter in the absolute magnitude, the 
scatter in the distance, and the photometric errors, it 
is likely that these stars have a relative scatter in their 
distance of the order of 25% or less. Again, this is a 
further hint that the "lower branch" stars might be a 
bulge population. 

— The distribution of "lower branch" stars in the color- 
color diagram (Fig. ^ forms an elongated region par- 
allel to reddening vector. This shows that at least part 
of the large scatter in their intrinsic colors might be 
due to reddening, local or along the line of sight. 

The items above suggest that the "lower branch" might 
be populated by evolved intermediate mass stars (about 
1 to 7 Mq), and likely on the Asymptotic Giant Branch 
(AGB), at about the distance to the Galactic center. We 
find a similarity between the "low er branch" stars and the 
sample of Galactic OH/IR stars in lJimenez-Esteban et alJ 
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Fig. 6. The color-color diagram for the selected stars in 
the field, after the removal of the set B2 of Eq. Q. A 
comparison with Fig. |21 shows that we were able to vir- 
tually remove all significant contamination from spurious 
reddening. 

l)2005|) . These authors conclude that the luminosity and 
color distribution observed could be explained as the result 
of a combined emission of a cool star (T ~ 2500 K) and a 
much cooler dust shell (T < 800 K). The main effect of the 
shell in the near-infrared colors would be to increase the 
circumstellar reddening, which, together with differential 
interstellar extinction effects from source to source, could 
explain the width of the distribution observed. Although 
we will not explore further in this paper the nature of these 
stars we note that they are particularly interesting objects, 
and given the large number of identifications in our fields 
(we have approximately 45 000 Bi stars and 1 500 B2 ones) 
a dedicated study is warranted. 

The "lower branch" stars seem to be unrelated to the 
molecular cloud, but their colors would be interpreted by 
the Nicer algorithm as a sign of extinction. This, clearly, 
could in principle bias our results toward a higher extinc- 
tion, especially at low galactic latitude regions. However, 
we argue that the bias introduced by "lower branch" stars 
is negligible. Indeed, the density of these objects is, in the 
worse case (set Bi) only 1% of the average density of stars, 
so that typically we have at most one "lower branch" star 
contaminating each pixel, and since we have ~ 30 stars 
per pixels, the effects of "lower branch" stars is negligible 
everywhere. 

Nevertheless, and in order to avoid any source of bias, 
although small, we excluded from the 2MASS catalogs all 
stars located in the B2 color-space region, and performed 
the whole analysis described in this paper using this re- 
duced subset of stars. We stress that if we had performed a 
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Fig. 10. The expected error versus the measured ex- 
tinction Ak for the various pixels in the extinction map 
of Fig. [71 shown as a density plot. 



a cut in the observed colors, we would be have introduced 
a new bias in the deduced column density; instead, the use 
of a select in the intrinsic colors does not bias the final 
results. As an example, Fig. shows the color-color dia- 
gram for the new set of stars: note that the "lower branch" 
disappears completely in this plot, a further confirmation 
that our selection is effective in removing this population 
of stars. 

We then run again the whole Nicer pipeline on the 
refined catalog. After (re)evaluating the statistical prop- 
erties of stars in the control field, we constructed the final 
map, shown in Fig. [71 in a grid of approximately 1 000 x 750 
points, with scale 30 arcsec per pixels, and with Gaussian 
smoothing characterized by FWHM = 1 arcmin; more- 
over, we used an iterative fi-clipping at i-a error. The 
final, effective density of stars is ~ 8 stars per pixel (this 
value changes significantly on the field with the galactic 
latitude, see Fig.lHJ; this guarantees an average (1-cr) error 
on Ax of only 0.019 magnitudes; the largest extinction is 
measured close to Barnard 59, where A^ — 2.68 mag (cor- 
responding to approximately Ay ~ 24 mag). As clearly 
shown by Fig. [71 the combination of the use of the 2MASS 
archive with the optimized Nicer technique allows us to 
reveal an unprecedented number of details. A quantitative 
analysis of this extinction map is delayed until Sect.0] 

Figure|51shows the expected error on Ak for each pixel 
of the extinction map. This figure deserves a few com- 
ments. First, we note that the most significant variations 
in the expected error are due to bright stars, which pro- 
duce the characteristic cross-shaped patterns (cf. Fig. [HJ . 
We observe that, with the help of Fig. [HI some apparent 
features in the extinction map are actually recognized to 
be due to the larger noise expected close to bright stars 
(these "features," hence, are just increased statistical devi- 
ations due to the increased local error) . This clearly shows 
that a detailed analysis of the extinction map of Fig. [3 is 
better carried out using also the density map of Fig.|Sland 
the error map of Fig. [51 Figure [HI shows also an increase 
in noise in the densest regions of the cloud, due to the 
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Fig. 7. The Nicer extinction map of the Pipe nebula. The resolution is FWHM = 1 arcmin; the contours are at 
Ak = {0.5, 1, 1.5} mags. 



reduced density of background stars. Finally, we note the 
increase of the error with the galactic latitude, also due to 
a change in the density of background stars. 

The accuracy of the column density measurements ob- 
tained in our field is shown by Fig. ^| which plots the 
expected error on Ak as a function of Ak for all pixels in 
our field. The exquisite data used allowed us to keep the 
average error well below 0.02 mag in Ak, and still to have 
a 1 arcmin resolution in our maps. The dynamical range 
of the Nicer extinction map can be better appreciated 
by noting that the lowest contour in Fig. [T] represents a 
10(7 detection, and that clumps such as Barnard 59 have 
a significance as large as ~ 200(t. Note that the increase 
on the error observed for Ak < 0.2 mag in Fig. ^| is due 
to the fact that regions with low extinction (including the 
control field) are located at high galactic latitudes, and 
hence have a smaller density of background stars. 

3. Distance 

An accurate determination of the distance of molecu- 
lar clouds is of vital importance to obtain a reliable 
estimate of the mass and of other physical properties. 



Unfortunately, distance measurements of molecular clouds 
are frequently plagued by very large uncertainties. A sim- 
ple method used often is based on the association between 
the cloud and other astronomical objects whose distance 
is well known. 'Onishi e t al.l (ll999^ associate the Pipe neb- 
ula with the Ophiuchus complex on the base of projected 
proximity and radial velocity and use the distance of the 
latter, (160 ± 20) pc, as the distance to the Pipe ijChinil 
.1981.1 . 

An alternative approach is based on the number counts 
of foreground sta rs. The method, used for example in 
lAlves et'all ()l998l) , exploits the large reddening produced 
by some clouds, which makes the identification of fore- 
ground stars relatively easy; then, galactic mo dels (e.g. 
iBahcall fc Soneiralll98d IWainscoat et alJll992|) are used 
to infer the expected number of stars (for each possi- 
ble cloud distance) inside the cone created by the cloud. 
Finally, the number of foreground stars observed is com- 
pared to the prediction and the distance of the cloud is 
inferred. Although this method is often the best one can 
use, it is unable to give accurate distances for several rea- 
sons: (i) it relies on galactic models, which might be in- 
accurate (especially at the small angular scales often used 
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Fig. 8. The star density map for the extinction map of Fig. \7\ i.e. the number of stars inside a (Gaussian) 1 arcmin 
beam. The dark structures correspond mostly to the position of dense cores in the Pipe nebula where the star density 
in minimal due to high extinction. The black "holes" outside the Pipe body correspond to masked bright stars while 
the white spots correspond to globular clusters. 



for molecular clouds); (ii) it is plagued by Poisson noise 
(because the number of stars inside the volume of the cone 
is a random variable) and (iii) stars do cluster (and thus 
the error is actually larger than the one expected from a 
Poisson statistics; imagine, in the extreme, the case of an 
unknown open cluster in front of the cloud). 

A more robust determination of the distance of 
the Pipe molecular complex can be obtained us ing the 
Hipparcos and Tycho catalogs ijPerrvman et alJ Il997^. 
The method, similarly to the star number counts described 
in the previous section, is based on the identification of 
foreground and background stars (observed on the line of 
sight of the cloud) for which a parallax estimate is avail- 
able. An upper limit for the distance of the cloud is thus 
given by the distance of the closest background star, i.e. 
the closest star showing a significant extinction in its col- 
ors. This novel approach to the distance of molecular cloud 
complexes has already been s ucces sfully applied to sev- 
eral clouds bv lKnude fc Hod l)l998l) . Here we revisit the 
method and use it to obtain a distance estimate for the 
Pipe nebula. 



We selected Hipparcos and Tycho-1 (1/239) stars ob- 
served in the area defined in Eq. Q, and matched these 
stars with the "All-sky Compil ed Catalog of 2.5 mil- 
lion stars" (ASCC-2.5, I/280A: iKharchenkd lioOl]) and 
with th e "Tycho-2 Spectra Type Catalog" (Tycho-2spec, 
HI/231: lmight et alJEool . The choice of these two lat- 
ter datasets (which are cross-references and merges of 
many stellar catalogs) was dictated by the need to ob- 
tain for each star with measured parallax an estimate of 
its spectral type. In particular, the use of ASCC-2.5 and of 
Tycho-2spec effectively allowed us to cover a large sample 
of spectroscopic catalogs: the Hipparcos catalog (1/239), 
the Carlsberg Meridian Catalogs (CMCll; 1/256), the 
Position and Proper Motions (PPM; 1/146, 1/193, 1/208), 
the Michigan Catalogs (III/31, III/51, HI/SO, III/133, 
III/214), the Catalog of Stellar Spectra Classified in the 
Morgan-Keenan System (III/18), the MK Classification 
Extension (III/78), and the FK5 catalog parts I and II 
(1/149 and 1/175). 

We considered all stars in our field with measured par- 
allax larger than the parallax error, and with spectral type 




Fig. 9. The map shows the statistical error (tai^ on the measured column density. Note that correlation on the errors 
on a scale of FWHM — 1 arcmin is expected. This map shows that for almost the whole field Ak — 0.05 represents a 
3-cr detection. 



B, A, F, G, K, or M. By comparing the expected B — V 
color (taken from Landolt-Bornstcin 1982, p. 15) with the 
observed one, we obtained an estimated of the color ex- 
cess E{B — V); we finally converted this into an extinc- 
tion in the V band A y by using a normal red dening law, 
Av = 3mE{B~V) JR.ieke fc T.ehofskJl flS.'j) . Note that, 
whenever possible, we used 2D spectral types; for ID spec- 
tral types we assumed a luminosity class IV. 

A plot of the star column density versus the Hipparcos 
and Tycho2 parallaxes is shown in Fig.^] Because of the 
relative large scatters in the parallax and column density 
measurements, Fig. llll is not straightforward to interpret. 
However, we note the following points 

1. At large parallaxes, for tt > 12 mas, the measured ex- 
tinction is low and consistent, within the errors, with 
a constant value. The average extinction in this range, 
{Ay) = 0.12 mag, either suggests the presence of a 
thin "veil" very close to us or of a small systematic 
error in the estimate of Ay. The presence of absorbing 
material towards the location of the Pipe nebula at 
distanc es below 50 pc is also suggested by the recent 
work of lLallement et al. (,,2003,1 . where they construct 



3D absorption maps of the local distribution of neu- 
tral gas through measurements of equivalent widths 
of the interstellar Nal D-line doublet towards nearby 
stars. In any case, from Fig. ^] we can deduce a ro- 
bust lower limit on the cloud distance, d > 80 pc. Note 
that, interestingly, we also detect a near single star 
(HD 158233) at (55.0 ± 1.67) mas with significant red- 
dening. Ay — (0.38 ±0.06) mag: thus, if the reddening 
has to be attributed to a thin molecular cloud, its dis- 
tance would be smaller than 20 pc. 

2. As we go to smaller parallaxes we start seeing stars 
with signs of a possible reddening. A first candidate 
is a star at 90 pc with Ay = (1.11 ± 0.66) mag (well 
visible in Fig. II 1(1 : however, the large error on Ay does 
not allow us to assess securely that the star has been 
significantly reddened by the cloud. 

3. At approximately tt — 7.3 mas (~ 137 pc) we ob- 
serve four stars with relatively large extinction. Ay > 
0.9 mag, and with relative small errors (the ratio 
Ay /Eri{Ay) ranges from 4.3 to 27.5). Hence, these 
stars give an upper limit to the cloud distance, d < 
140 pc. 
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Fig. 11. The reddening of Hipparcos and Tycho stars in 
the field. For this plot we selected only the stars charac- 
terized by a measured parallax it smaller than its error 
(Ttt, and have an estimated error on the column density 
S'Av < 1 mag. Among the f5f stars that matches these 
constraints, we marked with a filled dots the 44 ones that 
are located at 6 > 3°. 
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Fig. 12. The likelihood function of Eq. |7J as a function 
of the cloud parallax Trpipe marginalized over all the other 
parameters. The function shows a marked peak at tt ~ 
7.13 mas. The mean of Trpipc is 8.07 mas, while the median 
is 7.61 mas, corresponding to a distance of approximately 
130 pc. 



In summary, the simple argument discussed above sug- 
gests that the Pipe nebula is located between 80 and 
140 pc. 

A more rigorous estimate of the distance can be ob- 
tained by using a statistical model for the parallax- 
extinction relation. The physical picture here suggests 
that we can take the star column density as a random 



variable whose distribution only depends on the intrinsic 
measurement errors and on the parallax of the star. The 
simplest sensible approach is to consider a bimodal distri- 
bution for each star column density, where the extinction 
of stars in front of the nebula is consistent (within the mea- 
surement errors) with zero, and the extinction of stars at 
distances larger than the cloud can be either still zero (if 
the star is not observed through the cloud) or finite and 
with a large scatter (if the star is behind the cloud). In 
this case, the conditional probability of observing a star 
with visual extinction Ay given the fact that the star is 
located at a parallax tt is 

r Gau(Ay I A^f , af^ + if tt > ^Pipe 

[ + /Gau(Ay|<s, aX + if ^ < ^Pipe 

(6) 

where we denoted with Gau(a;|a;, cr^:) the value at x of 
a normal probability density with mean x and variance 
(T^ [cf. Eq. H12|) ] . In other words, we use a normal dis- 
tribution with low average A^y for foreground stars, and 
the mixture of the same distribution with another normal 
(with large average A^) for background stars. For both 
distributions the variance was set to the sum of <j\^ , the 
estimated variance in the measured column density of the 
star, and crj^ (respectively, c^^), the intrinsic variance 
of the foreground (background) distributions. The param- 
eter / G [0, 1] used in Eq. © represents the filling factor, 
i.e. the ratio between the area of the sky occupied by the 
cloud and the area of the whole field. 

We considered this simple model leaving as free param- 
eters {TTpipe, /, Al^, } (the average foreground Ay and 
its variance were deduced from the stars with measured 
parallax tt > 12 mas). In order to assess the goodness of 
a model, we computed the likelihood function, defined as 



TTlplTT TT 



) dTT 



(7) 



where the product is carried over all N observed stars with 
parallax, and where pi^n | tt^"-*) is the probability distri- 
bution that the n-th star with measured parallax -n"*^"^ has 
a true parallax tt. The integral in Eq. ((TJ takes into ac- 
count the possibility that a star with a measured parallax 
< Ti'Pipc (or TT > TTpipc) is iucorrcctly taken as back- 
ground (respectively, foreground) because of its large par- 
allax error aj^. In our analysis, we used for p(^Tr | tt'^"^) a 
simple normal distribution with standard deviation equal 
to the formal error on the measured parallax ctI"^ : 



P[tt tt 



Gau(7r TT 



(8) 



We note that this particular choice allows us to evaluate 
the integrals of Eq. iQ analytically in terms of the error 
function erf. 

In order to study in detail the likelihood function (Q 
in its multidimensional parameter space we used Monte 
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Fig. 13. The distribution of individual star extinctions. 
The solid histogram shows the distribution of column den- 
sities for stars observed in low-extinction regions (—0.05 < 
Ak < 0.05 on the extinction map); the dashed line 
the same distribution for stars in high-extinction regions 
(1.15 <Ak< 1.25). 
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Fig. 14. The dispersion on the extinction measurements 
as a function of the column density Ak- In order to have 
a clear plot, we evaluated the dispersion on bins of 0.05 
magnitudes in Ak; the plot shows the average values ob- 
tained in each bin and the relative scatter. The dashed 
line represents the expected increase in the dispersion if 
we attribute this effect to the contamination of foreground 
stars. 



Carlo Markov Chains (MCMC; see, e.g. iTanned Il99lt) . 
The obtained results are summarized by Fig. I12L where 
we show the marginalized likelihood as a function of the 
cloud parallax tt. As best estimate for the cloud distance 
we take the median of the distribution shown in Fig. ^1 
rfpipo = 130 pc; the formal 68% (respectively, 95%) confi- 
dence regions are [110, 143] pc and [72, 154] pc. Note that 
the best fit value obtained for the filling factor is / = 0.42, 
which compares well with the value directly measured in 
the map of Fig.[7|(41% of the area in our field has a mea- 
sured extinction Ay > 3 mag). 

In summary, in the rest of this paper we will use a es- 
timate of the cloud distance 130^20 P*^- '^^^ formal error 
is relatively small, but it is probably underestimated (e.g., 
the error due to a possible misclassification of the spectral 
type of a star was not included in the error budget) ; more- 
over, because of the relatively small number of reddened 
stars, the distance obtained appears to be slightly model- 
dependent. Our distance estimate is likely to be biased 
towards large values because the method used gives an 
upper limit (this is also suggested by the long tail at large 
parallaxes in the likelihood function of Fig. ll2|) . With the 
advent o f the new generation astrom etric missions such 
as Gaia l)Lindegren fc Perrvmanl ll996') it will be possible 
with this or similar methods to accurately measure the 
distances of a large sample of molecular cloud complexes. 

4. Statistical analysis 

4.1. Star extinction 

In Fig. El we plot the distributions of column densities ob- 
tained for stars observed in low-extinction (—0.05 mag < 
Ak < 0.05 mag) and high-extinction (1.15 mag < Ak < 



1.25 mag) regions. The two histograms show two simi- 
lar Gaussian shapes, but the one corresponding to high- 
extinction has a larger width (the best-fit Gaussian dis- 
persions are cti — 0.210 mag and 172 — 0.272 mag). This 
increase can be due to a number of factors: (i) the in- 
crease of photometric errors for the faint stars observed 
through dense clouds, and the lack of the K band photom- 
etry for most of these objects; (ii) the internal structure 
of the cloud on scales smaller than our resolution (which, 
we recall, is FWHM = 1 arcmin); (iii) an increase in the 
contamination of foreground stars (see below Sect. 14. 41) . 
Among these factors, (i) can be evaluated from the pho- 
tometric errors of the 2MASS database, and appears to 
negligible (i.e., there is no significant increase in the pho- 
tometric error of stars in moderately extincted regions, cf. 
Fig. I10|l . On the other hand, factors (ii) and (iii) are dif- 
ficult to disentangle, because they essentially produce the 
same effects (see below). 

The increase of the dispersion on Ak at large col- 
umn densities can be evaluated more quantitatively from 
Fig. El This figure has been obtained from the map of 
the dispersion of Ak (not shown here). This map is con- 
structed by evaluating, for each pixel of Fig. |3 the scat- 
ter of the column density estimates corresponding to each 
star. Hence, large values in some pixels of this map mean 
that the stars used to estimate Ak in the correspond- 
ing pixels of Fig. |7| present very different reddening val- 
ues (see Paper I for the analytic definition of the scatter 
map). Figure ITU has been obtained by plotting the aver- 
age values observed in the dispersion map in pixels char- 
acterized by a given value of Ak in Fig. [7| We note that 
the observed increase in the dispersion of Ak has been 
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observed in several similar studies (e.g. iLada et al.lll994t 
lAlves et al.lll998l iLada et al.lll999(l . 

As we mentioned above, in our case it is dilficult to 
attribute the effect observed in Fig. with certainty to 
foreground stars or to substructures. Consider, as an ex- 
ample, a "two-dimensional Swiss cheese" molecular cloud, 
i.e. a thin cloud with many holes (some of them smaller 
than the resolution of our maps). In this case, a back- 
ground star observed through a small hole will produce 
exactly the same statistics as a foreground star. Hence, in 
this case it will be impossible to distinguish foreground 
stars and substructure. In less extreme (and more phys- 
ical) situations it is generally possible to firmly identify, 
in the dense regions of the cloud complexes, foreground 
stars. However, this task is unfortunately non-trivial for 
the Pipe, because of its filamentary structure and of the 
relatively low resolution attainable with the 2MASS data. 
Still, our estimate indicates that only an extremely small 
fraction Fq ~ 0.003 of foreground stars is present in our 
map (see below Sect. l4.4|l . and thus they have a negligible 
contribution to the effect observed in Fig. (cf. dashed 
line in that figure) . In summary, we attribute most of the 
increase in the dispersion of Ak to unresolved substruc- 
tures. A further indication of this is given by similar anal- 
yses of molecular clou ds carried out at a higher resolution 
(e.g. iLada et alJEoO^ , which show a dispersion sub- 
stantially independent of Ak (i.e., a flat curve in the Ak- 
(JAk plot). In a follow-up paper (Lombardi 2005hl we will 
further investigate the effect of substructures in molecular 
clouds. 

4.2. Reddening law 

The reddening law plays a fundamental role in any in- 
frared dust measurement since it is used to derive the 
column dust density from the properties of background 
stars. In particular, a key assumption is that the redden- 
ing law is linear, i.e. that the ratio A\.^/Ax^ of extinctions 
in two different wavelengths is constant. Since the intrin- 
sic luminosities of stars are not known, we cannot measure 
directly the extinction Ax in a single band; in the NIR, 
however, we can measure with good accuracy the color 
excess E{m\-^ — mx^) = Ax^ — Ax2 , and this is the essence 
of the Nice and Nicer techniques. Moreover, if we have 
good photometry in three bands, say J, H, and K, we 
can indirectly verify the assumption of a linear redden- 
ing law by checking that the color excesses E{H — K) 
and E{J — H) are linearly dependent. The large number 
of stars in the field allows us to perform this check with 
great accuracy. 

The linearity of the reddening law, ultimately, is equiv- 
alent to say that the stars in the color-color plot are mostly 
found along a linear stripe, and the hypothesis is qualita- 
tively well supported by Fig. IHl Still, this simple observa- 
tion is not easily translated into a robust and statistically 
accurate method to measure the reddening law. This, in 
part, might be the origin of some discrepancies found in 
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Fig. 15. The reddening law as measured on the analyzed 
region. The plot shows the color excess on J — if as a 
function of the color excess on H — K (the intrinsic colors, 
deduced from the control field, are 0.18 and 0.60 respec- 
tively); the solid line shows the best fit. Error bars rep- 
resents the standard deviation of color excesses inside the 
bin. 



the literature about the slope of the reddening law (an- 
other good reason are the different NIR photometric sys- 
tems used by different authors). 

For example, in order to measure the slope of the red- 
dening law, we can fit all star colors (or color excesses) 
with a linear relation: 

[J -H) = a{H -K) + f3 , (9) 

where a — {Aj — Ah) /{Ah — Ak) and (3 are taken to 
be constant for all stars. Naively, the fit could be done 
directly on the individual star colors {Hn — Kn} and {Jn — 
Hn}i by minimizing the expression 

N 

X^^Y. ["(^" ~Kn)+f3- {Jn - Hn)] ' , (10) 
n=l 

or similar expressions taking into account the photometric 
errors of stars. In reality, such a procedure leads to unsat- 
isfactory results because the large majority of stars are 
subject to small or negligible reddening. In other words, a 
minimization of Eq. (|10|l would basically take into account 
only the lower-left part of the plot of Fig. El and would 
neglect the most interesting (highly reddened) stars (say, 
H - K > 1 mag). 

In order to equally take into account all reddening 
regimes, we followed here a different procedure: 

1. We evaluated the column density of each star by us- 
ing the normal reddenin g law considered in this paper 
l|Tndebetouw et aljEool . 

2. We binned the stars according to their extinction in 
bins of 0.02 mag. 
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3. In each bin, we evaluated the average color excesses 
in H — K and J — H, and the standard deviation in 
J — H . Note that the standard deviation in J — H was 
calculated by referring all star colors inside the bin to 
the average H — K color, or equivalently by evaluating 
the standard deviation of the expression 

E{J - HUrr = - a[E{H -K)- {E{H - K))] 

+ E{J-H). (11) 

This procedure ensures that the scatter in E{J—H)coii 
does not include the scatter due to the finite size of the 
bin inK-H. 

4. We fitted all positive bins with a linear relationship of 
the form of Eq. (in the fit we took into account the 
standard deviations evaluated in the previous point). 

5. Finally, we repeated the whole procedure using the 
newly determined reddening law until convergence is 
achieved (which typically happens after two or three 
iterations). 

The method presented above has many pleasant proper- 
ties. First, all reddening regimes equally contribute to the 
fit, since by construction we equally use all bins of different 
column densities. This, in turn, implies that the method 
is robust against bright stellar populations that might ar- 
tificially change the slope of the reddening law (cf. above 
Sect.|2). Moreover, the use of the observed standard de- 
viations on the J — H color for the fit properly takes into 
account local deviations from a reddening law or increased 
photometric errors for reddened stars. 

We applied the method described above to our data, 
by selecting all 2MASS stars with accurate photometry in 
all bands (we required all photometric errors to be smaller 
than 0.1 mag) from the cleaned catalog (cf. above Sect. [21); 
the results obtained are summarized in Fig. 1151 We found 
as best fit slope a = 1.82 ± 0.03, a v alue that appears 
to be in excellent agreement both the llndebetouw et all 
\2m^ 2MASS reddenin g law, E{J - K)/E(H - K) = 
1.78 ±0.15, and with the lRieke fc Lebofskvl l)l985|) E{J - 
K)/E{H — K) — 1.70 normal reddening law converted 
into the 2MASS internal photometric system (iCarpenteil 
l200i- which is E[J - K)2UAss/EiH - if)2MAss - 1.89. 

In follow-up papers we will apply the same technique 
to other cloud complexes studied from the 2MASS archive. 
The uniformity of the 2MASS data and of the procedure 
used to derive the reddening law will allow u s to accurately 
study cloud-to-cloud variations (see, e.g., iKenvon et al.1 
l|l998l) for a case where significant differences from a stan- 
dard reddening law are found). 

4.3. Smoothing method 

As discussed in Paper I, the pipeline developed allows 
us to obtain smooth maps of the discrete column den- 
sity measurements (obtained for each background star) in 
three different ways: simple average, sigma-clipping, and 
median map. For the simple average we use all column 
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Fig. 16. Plot of the column densities for the various pix- 
els of the extinction map obtained through median fil- 
ter and simple average. The diagonal line shows the locus 
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Fig. 17. This plot, similarly to Fig. 1161 shows the rela- 
tionship between the pixel column densities measured us- 
ing the median and the simple average. The plot has been 
produced by splitting the pixels in bins of 0.05 magnitudes 
in ^median ^ ^.^^ plotting, for each bin, the average and 
the 2-CT scatter of the quantity ^^''^'^ - yl^odian^ 



density measurements close to each pixel, and the extinc- 
tion value assigned to the pixel is a weighted average of 
the Ak measurements for each star. The weight takes into 
account both the statistical error on Ak (which is due to 
the photometric errors and to the intrinsic scatter of star 
colors), and the angular distance between the star and 
the pixel. For the sigma-clipping, instead, we selectively 
discard outliers on Ak, i.e. stars that have associated col- 
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Fig. 18. This figure compares the column densities ob- 
tained through sigma-cHpping and simple average, simi- 
larly to what done in Fig. El for the median. 



umn densities significantly different from the local aver- 
age. Finally, the median map is constructed by evaluating 
a sort of "weighted median" for the column densities mea- 
sured for angularly close stars. 

Figure ll()l shows the relationship between the extinc- 
tion measurements obtained through simple average and 
median. The very tight band show that the two methods 
give comparable results. A more quantitative analysis is 
provided by Fig. El where an histogram of the differences 
between the simple mean and the median for each pixel 
is presented. Note that the two methods are statistically 
indistinguishable up to Ak — 1.5 mag, where the mean 
starts to estimate consistently smaller column densities. 
Hence, this confirms that the use of the simple mean is jus- 
tified up to relatively high column densities. FigurelTTIalso 
shows an increase in the scatter at high column densities, 
an indication of a difference in the extinction estimates 
at the corresponding sky locations using the mean and 
the median smoothing techniques. Unresolved structures 
in the cloud, which are expected in the dense regions, are 
the most probably reason for the observed scatter (sub- 
structures are likely to play different roles in the simple 
mean and median smoothing). 

A statistical comparison between the simple mean and 
the cr-clipping smoothing techniques is shown in Fig. El 
From this figure, we can deduce that there is no significant 
difference in using these two techniques. This is mostly due 
to the fact that the two methods are expected to provide 
different results only at high Ak, when a possibly signif- 
icant number of foreground stars might contaminate the 
map. On the other hand, these regions also have a much 
smaller star density, and thus the extinction map there 
will have a large error (see Fig. jSJ. This, in turn, implies 
that the sigma clipping is not effective there (the statis- 
tical error on the map is large enough to accommodate, 
within a 3-cr interval, unreddened stars). Hence, in the 
following we will only consider the median and the simple 
mean estimator. 
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Fig. 19. A zoom of the extinction map showing 
Barnard 59. Note the white "hole" close to the center of 
this cloud, which is due to the combined effect of very high 
extinction values and the presence of embedded stars in 
this star forming core. 



In conclusion, our analysis shows that the simple mean 
estimator, which has the smallest scatter, is reliable up to 
approximately Ak ~ 1.5 mag; for larger column densities, 
the median should be used in order to minimize the effect 
of foreground stars. 

4.4. Foreground star contamination 

The three smoothing methods described above differ sig- 
nificantly on the effect of foreground contamination. If we 
assume that a fraction F > of stars arc foreground with 
respect to the cloud, then the simple average estimate will 
be biased toward small column densities. In particular, the 
average measurement of Ak will be (Ak) = (1 — F)A'^^^°. 
In contrast, the median will provide an ( almost) unbiased 
estimate oi Ak as long as F < 0.5 (see ICambresv et all 
120021 and lLombardill2005al for a detailed discussion of the 
median properties). The sigma-clipping estimate will of- 
ten be between these two extremes: it will be effective in 
removing foreground stars only in relatively dense regions 
and only for small values of F. 

Because of selection effects, the value of F changes 
significantly on the field, and in particular increases in 
high-column density regions. This, in turn, implies that 
the correction to be used on the estimated value of Ak is 
not constant on the field. In the case of the Pipe nebula 
only a few magnitude extinction in V is observed in most 
of the field (see above Sect. 12 see also Fig. (23 below), 
and thus the value of F is expected to be approximately 
constant. Moreover, since this cloud complex is very close 
to us and is observed close to the galactic center, we expect 
only a tiny fraction of foreground stars. 
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Fig. 20. The distribution of pixel extinction over the 
whole field. The figure shows the histogram of column 
densities for all pixels (gray line) and for the pixels with 
galactic latitude & > 3°. Note that the latter histogram 
can be fitted almost perfectly using four normal distri- 
butions (solid line) or four lognormal distributions (not 
shown here), while a single lognormal (dashed lined) is 
able to reproduce the general shape of the histogram. 



In order to evaluate quantitatively the density of fore- 
ground stars, we have selected high-extinction regions 
characterized by Ak > 0.8 mag and with relatively small 
expected error in Ak (we allowed for a maximum error of 
0.2 mag). In order to avoid the effects of possible substruc- 
tures and ambiguities in the identification of foreground 
stars, we restricted our analysis to the main core of the 
Pipe nebula, i.e. to the region at I ~ 1.5° and 6 ~ 4° (see 
Fig. d). We have then checked all stars in these regions 
that show "anomalous" extinction, i.e. stars whose col- 
umn densities differ by more than 3-(t with respect to the 
field. A total of 70 stars met this criteria; hence, since the 
area selected is approximately 1057 arcmin^, we estimated 
a foreground star density of 0.066 arcmin^^. As a result, 
the relative fraction of foreground stars in regions with 
negligible absorption is only F\av=o = Fq = 0.003, and 
we can safely ignore the effect of foreground stars except 
on the higher extinction regions. Note, in particular, that 
a clear sign of very high extinction (very few background 
stars) and contamination by foreground stars is observed 
only in Barnard 59, where we see a "hole" in the extinction 
close to the center of this clump (see Fig. ll9|) . In this case, 
these foreground stars are not really foreground but young 
stars moderately em bedded in Barnard 59, an active star 
forming region (e.g. lOflishrii^lliii)- 

Our estimate of 

Fq allows us to evaluate the maximum theoretical extinc- 
tion meas urable with the Nic e and Nicer method: using 
Eq. (53) of lLombardil l|2005a(l . we obtain yl^^'^ ~ 2.9 mag, 
which is close to the maximum value measured here. 

4.5. Column density distribution 

A fundamental statistical pro perty of a cloud is the dis- 
tribution of column densities. IVazouez-Semadenil l)l994|) 
showed that for highly supersonic flows (such that the 



Mach number AI — u/cs ^ 1), the gas essentially has 
a pressureless behavior, and gravitational forces can also 
become negligible. In these conditions, the hydrodynamic 
equations become scale- invariant: in other words, mo- 
tions at different length and density scales obey essen- 
tially the same equations. For a fully developed turbu- 
lent gas, density and velocity can be regarded as random 
fields. Because of the scale invariance, the probability of 
having a local (volume) density fluctuation of amplitude 
Ap depends uniquely on the ratio Ap/p, i.e. on the rel- 
ative fluctuation amplitude; hence, the probability distri- 
bution function of the density at each point of the cloud 
is expected to be lognormal. When projecting the three- 
dimensional mass density along the line of sight, the log- 
normal distribution is essentially preserved, i.e. the pro- 
jected two-dimensional density is also expected to be well 
approximated by a lognormal distribu tion. 

More recently, IVazauez-Semadenil argument has been 
investigated in more detail u sing various simulatio ns 
by a number of authors fe.g. IPadoan et all ll997Jbl lal: 
IPassot fc Vazauez-Semadenilll998HKlessenll200nl) . and in 
all cases the lognormal distribution was found to be a good 
approximat i on of the (projected) cloud density. However, 
IScalo et all l)l998j) found that, contrary to previous claims, 
the density probability distribution is well approximated 
in their simulations by an exponential law in a relatively 
large range of physical parameters. Re cently, this argu- 
ment has been challenged by Ostrikc r et alJ ||2Q011 by 
showing a good agreement wit h a lognormal dis tribution 
in the molecular cloud IC5146 l)Lada et al.lll999(l . 

A direct study of the extinction distribution on a real 
cloud is a stringent test for simulations and ultimately also 
provides important hints on the physical conditions inside 
molecular clouds. We evaluated the column density distri- 
bution in the Pipe nebula by constructing the histogram of 
extinction measurements in each pixels of the map shown 
in Fig. 13 with 6 > 3°, binned on suitable intervals (we 
used here 0.002 mag in Ak). The obtained distribution, 
plotted in Fig. 1201 shows a complex shape with multiple 
peaks. 

In order to compare the empirical distribution with a 
theoretical one, we used Poisson models for the Individ- 



Distr. 






o 


rel. factor 




(mag) 


(mag) 


(mag) 




InGau #1 


-0.241 


0.503 


0.325 


1.000 


InGau #1 


-0.150 


0.413 


0.501 


0.150 


InGau #2 


-1.805 


1.943 


0.018 


0.744 


InGau #3 


0.039 


0.267 


0.371 


0.106 


Gau #1 


0.015 




0.022 


0.086 


Gau #2 


0.123 




0.033 


0.353 


Gau #3 


0.267 




0.064 


0.395 


Gau #4 


0.419 




0.118 


0.166 



Table 1. The best-fit parameters obtained from a fit of the 
column density distribution shown in Fig. l20l using a single 
lognormal distribution, three lognormal distributions, or 
four normal distributions. 
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ual histogram bins, and evaluated the joint probability 
to have the observed histogram in the range —0.3 mag < 
Ay < 2 mag (we excluded high Ax bins in order to avoid 
the complications inherent with large-column density re- 
gions). We constructed the theoretical models as sum of 
normal distributions 



Gau(a;|a::o, cr^) 



1 



V27 



: exp 



{x - XoY" 
2cr2 



(12) 



and lognormal ones 

lnGau(a;|xo, xi, cr) = 
1 



: CXp 



[ln(a: — xq) — lna;i] 
2^ 



(13) 



As shown by Fig. 1201 the best fit obtained from a sin- 
gle lognormal (dashed line) is surprisingly good given the 
complexity of the distribution; clearly, however, in this 
case we are only able to reproduce the general shape of 
the histogram and we are unable to generate finer details 
clearly visible in Fig. 1211 As a measure of the quality of the 
fit, we note that the residuals are more than ISOcr away 
from the observed histogram. 

In order to improve the fit we added more components. 
However, two-component distributions still gave unsatis- 
factory results, and we had to use more components distri- 
butions to obtain a close match to the data. In particular, 
we verified that both the sum of four normal and three 
lognormal distributions gave similar fits, with residuals of 
the order of 8a. Hence although, in principle, the fit is 
still statistically inconsistent with the data, we note that 
in practice the large number of bins and the non-trivial 
structure of the Pipe nebula make it difficult to obtain 
good fits with a relatively small number of parameters. 
Note also that the need for several components in the 
column-density distribution is reminiscent of the veloc- 
ity structures obser ved in the Pipe nebula from CO data 
llOnishi et a,1.lll99!ll) . In this respect, it is also likely that 
some of the "weak" components (e.g., the InGau #3 or 
the Gau #4) are related to background clouds observed 
in projection to the Pipe nebula at low galactic latitudes. 

5. Comparison with CO analysis 

Radio observations of H2 surrogates, and in particular of 
CO isotopes, provide an alterna tive independent est imate 
of the cloud gas column density. lOnishi et alJ ( 199fll) stud- 
ied in detail the Pipe nebula with the NANTEN radiote- 
lescope and kindly provided their ^^CO map to perform a 
comparison with the Nicer analysis described in this pa- 
per. To this purpose, we constructed an extinction map us- 
ing the same resolution (4 arcmin) and coordinate system 
of the NANTEN observations. In order to avoid any con- 
tamination by low-galactic latitude clouds, we excluded 
all measurements at 6 < 3° (see Fig. 1^. 
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Fig. 22. The -^^CO integrated intensity /i2co as a function 
of the Nicer extinction. Each p oint of this plot show s the 
^2 CO radio measurement from lOnishi et al and 
the Nicer extinction at the corresponding position. The 
filled squares represent the average ^^CO intensity in bins 
of 0.1 mag of X-band extinction; the solid, smooth line 
is the best fit for functions of the form of Eg. II14II; the 
dashed line shows the linear relation used bv lOnishi et alJ 



(>199ai (X-factor 2.8 x 10^ 



K"^ km-i s). 




Fig. 23. The residuals in the best-fit ^ of Fig.ES The 
filled squares represent the standard deviations in bins of 
0.1 mag. 



Figure [221 shows the relationship between the Nicer 
column density and the NANTEN integrated ^^CO tem- 
perature, /i2co = J Tdv. The large number of indepen- 
dent measurements shown in this figure (~ 5 000) was used 
to compare these two estimates of the gas column density 
and to show the limits and merits of both. 

From Fig. 1221 we can deduce a number of qualita- 
tive points. We first note that apparently ^■^CO mea- 
surements are insensitive to low column density regions. 
Specifically, it appears that up to ^ 0.2 mag of K-hand ex- 
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Fig. 21. Left: The NANTEN integrated ^^CO column density map (kindly provided bv lOnishi et alJll999l) : the white 
regions have not been observed and no data are thus available there; the shaded region is located at 6 < 3° and has 
been excluded from the analysis to avoid contamination from low-galactic latitude clouds. Right: The Nicer extinction 
map downgraded to the resolution of the NANTEN map; shaded regions are excluded from the analysis. 



tinction the radio measurements are uniformly distributed 
around K km s~^. The relatively small scatter observed 
in Fig. 1221 up to 0.15-0.2 mag of extinction (cf. also 
Fig. 123) indicates that both methods considered here have 
small intrinsic internal errors. 

At higher extinctions, and up to Ak — 0.6 mag, we 
observe an almost linear relationship between the extinc- 
tion and the ^^CO measurements; correspondingly, the 
scatter in the plot increases significantly (see below). We 
note that the presence of a linear relationship between 
the CO integrated temperature and the NIR extinction 
(and thus the hydrogen projected density), although on 
a relatively small extinction range, is not obvious on the 
scal es considered here and for non-virialized cloud systems 
Csee lDickman et allllQSi) . Finally, for > 0.7 mag, the 
"'^^CO data appear to saturate to a constant value close 
to 30 K km s^"'^. This well-known saturation effect is de- 
scribed in terms of an exponential relation between the 
integrated temperature h-^co oc 1 — e""^ and the cloud 
optical depth r cx A^- The following analysis will mostly 
focus on the unsaturated ^^CO regime and will thus di- 
rectly make use of the integrated temperature. 

A more quantitative analysis of Fig. |211was carried out 
as follows. We divided the measurements in regular bins 
in Ak (we used a bin size of 0.1 mag), and we computed 
in each bin the average of the ^^CO intensity. The results 
obtained are shown as filled squares in Fig. [221 This sim- 
ple plot confirmed the qualitative remarks discussed above 
and suggested that we could approximate the Ay-^^CO 
relationship with a function of the form 



2CO 



A 



1 



1 + cxp -(Ar- 



(14) 



We fitted this equation to the data by minimizing the 
scatters between the predicted CO integrated intensity 



2D 




Fig. 24. The ^^CO-Ak relation, with datapoints binned 
along the CO axis every 5 K km s~^ (filled squares). The 
solid line represents the best fit (on the whole field) from 
Eq. fT5|l . 



and the observed one; the best fit parameters obtained 
were A = 32.3 K km s"\ Af'^ = 0.51 mag, k = 
6.20 K km s~^ mag~^, and h — 0.036. The residuals of 
this fit with the data are shown in details in Fig. [23 
the increase of the dispersion in the relation (|14|l at 
Ak — 0.2 mag is evident from this plot. Since the ex- 
pected error in the Nicer map of Fig. [21 is as low as 
~ 0.01 mag, and since the expected error in the ^^CO 
integrated velocities is also relatively small (this can be 
estimated from the residuals at Ak — mag of Fig. 1231 
and is of order of 1.5 K km s~^), we can deduce that the 
scatter shown in Fig.[2Slfor Ak > 0.2 mag is physical: the 
ratio of dust and ^^CO in the Pipe (and likely in other 
molecular clouds) is far from being constant. 
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10 20 

Fig. 25. The scatter of our measurements on the linear 
fit of Eq. (|15|l . Note how the derived standard deviation 
(filled squares) is practically constant over the CO column 
densities investigated here. 

So far we investigated the Ak-^^ GO relationship using 
the value of Ak as independent quantity: in other words, 
we studied the expected CO measurement for each given 
Ak column density. We now swap the role of Ak and 
CO, and consider the average Ak value corresponding to 
a given ^^CO measurement. To this purpose, we averaged 
the values of the Nicer extinction in bins of 5 K km . 
The result, shown in Fig. suggests that we can well 
approximate the average with a linear relationship of the 
form 



A 



K 



A 



(0) 
K 



rll: 



CO 



(15) 



Note that we need to include explicitly a non-vanishing 
"zero point", A^^\ for the Ak measurement. This is due 
the dissociation of the CO molecule by the interstellar UV 
radiation field. Our results indicate that CO molecules in 
the Pipe become (self-) shielded from the interstellar ra- 
diation field at about 1 magnitude of visual extinction 
(2 magnitudes along the entire line of sight through the 
cloud), consistent with sta ndard theoretical predictions 
and prior observations (e.g.. Ivan Dishoeck fc BlacklllQSal 
lAlves et al.lll99<i iBergin et alJl2n02l) . This CO threshold 
should in principle be a function of the intensity of the lo- 
cal interstellar radiation field and could in principle vary 
from cloud to cloud. We stress that it is highly unlikely 
that the Nicer technique overestimates the extinction at 
low Ak, i.e. that the "zero point" observed in the rela- 
tion (|15|l is an artifact; rather, if there is a bias in Nicer, 
this is likely to be toward an underestimate of the column 
density (because of a possible reddened control field). 

Our data, combined with the ^^CO data, allows for 
the best determination of the CO-to-H2 conversion fac- 
tor (X- factor) using dust as a tracer of H2, because of 
the large number of measurements (approximately 5 000) 
and also because our (smoothed) dust extinction measure- 
ments have a mean error smaller than 0.05 magnitudes of 
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Fig. 26. The difference between the two maps of Fig. 1211 
i.e. the Nicer extinction map, converted into intensity 
using Eq. ((TH) . and the ^^CO intensity. 



visual extinction. To derive the X-factor for the ^^CO data 
we performed several best fits using Eq. l(T^ using differ- 
ent selections of the points of Fig. |^ the results obtained 
are reported in Table [3 

Finally, Fig. 1251 shows the residuals obtained from the 
fit of Eq. l(TH|l . Interestingly, these residuals appear to be 
rather constant on the field, and their standard deviation 
is approximately 0.1 mag in iiT-band extinction (or about 
1 magnitude of visual extinction): hence, a single pointing 
in "'^^CO can only constrain the dust column density with 
an accuracy of a few magnitudes. Remarkably, this re- 
sult is consistent with that found bv lLada et a l.' ("l994'l for 
"CO in the molecular cloud IC 5146.lAlyes et al.. tl99il 
for i^CO in L977, and lBergin et all l|2n02l) for the molec- 
ular cloud Barnard 68. Hence, we find what it seems to be 
an irreducible uncertainty in the dust-CO correlation of 
about 1 magnitude of visual extinction. The cause for this 
irreducible uncertainty is not known and it deserves a ded- 
icated study. A possible origin for this uncertainty might 
lie in the sensitivity of the CO molecule (and its isotopes) 
to variations in temperature (even along t he same line-of- 
sight ) , density and UV radiation fi eld fe.g. lFrerking et all 
ll982Hvan Dishoeck fc Blacklll988(l . We note however that 
changes in the abundance of ^^CO, due to freezing out on 
grains, are probably not in place at the extinctions less 
than Ak — 1 mag considered here when comparing the 
NANTEN and 2MASS/NiCER analyses. 

In order to further test this point, we considered the 
difference between the Nicer extinction and the ^^CO 
integrated intensity. In particular, we converted the de- 
graded Ak Nicer extinction map shown in Fig. |^ (right) 



into a ^^CO intensity using Eq. and subtracted from 
this map the NANTEN ^^CO extinction map. The result 
obtained, shown in Fig. 1261 gives insight into the origin 
of the scatter observed in Fig. The first striking fact 
is that clear structures are observed in Fig. 1261 and this 
alone rules out the possibility that the difference between 
the 2MASS-N1CER and the NANTEN measurements is 
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Fit A'^^ 

(mag) (mag km~^ s) 
Whole sample 0.176 0.0184 
/i2co > 5 K km 0.199 0.0172 
Ak > 0.1 mag 0.211 0.0165 
e [0.1,0.6] mag 0.226 0.0127 

Table 2. Best fit parameters relative to Eq. H15|) . 
If a normal reddening law is assumed, the X-factors 
derived from these fits are {4.21,3.93,3.77,2.91} x 
10^° cm^^ K^^ km^^ s respectively. The formal, la er- 
rors are 0.002 mag and 0.0002 for A^' and r, respectively 
(for all four fits). 



due to statistical errors. Instead, in large areas we observe 
systematically positive or negative values of the difference 
between the two maps. 

In analysing Fig. 1261 one should always keep in mind 
that this map has been built by converting the Nicer 
column density into a radio intensity by using the best fit 
provided by Eq. H14I) : hence, differences have to be inter- 
preted as deviations from the fit used. With this point in 
mind, we can deduce some interesting facts from Fig. 1261 

— The outer regions of the nebula appear to have consis- 
tent measurements in radio and NIR. Given the form 
of Eq. (|14|1 . this further confirms that ^^CO observa- 
tions are insensitive to low column densities. 

— In many peripheral parts of the Pipe nebula, the 
^^CO intensity appears to overestimate the (con- 
verted) Nicer column density (see, e.g., the dark re- 
gion in the middle of Fig. [211) • This effect, in principle, 
could be related to differences in the photo-dissociation 
of CO molecules due to variations of the local interstel- 
lar radiations field. However, since generally an excess 
of ^^CO intensity is observed in the periphery of the 
Pipe, it is more likely that the effect is related to tem- 
perature gradients in the nebula (this is also suggested 
by the patchy aspect of the map of Fig. ES|l . 

— The dense cores (e.g., Barnard 59) generally appear 
bright in Fig. |2S1 indicating that the ^^CO map un- 
derestimates the column density there when compared 
to the Nicer analysis. This effect is clearly due to the 
saturation effect of ^^CO at high column densities (the 
saturation hence appears to be even stronger than sug- 
gested by Fig. 121 and by the fit ifTH) ). 

The items above proves that a joint analysis of the Nicer 
and CO results can provide nontrivial and interesting in- 
sights into the structure of molecular clouds and the pro- 
cess that regulate the radio emission. Clearly, the results 
obtained here are specific to the Pipe nebula; however, by 
carrying out similar discussions for other cloud complexes 
we hope to be able to draw general conclusions on this 
important topic. 



6. Mass estimate 

The cloud mass M can be derived from the Ak extinction 
map using the following simple relation 

M = (ffi(3K [ Ak (Px , (16) 
Jn 

where d — 130 pc is the cloud distance, fi is the mean 
molecular weight corrected for the helium abundance, 
13k ^ 1.67 X 10^^ cm~^ mag"^ is the rat io jV(Hi) + 
N(li^.)/AK JSavage fc Mathijll979t see also EhU 
iBohlin et al.l Il973) . and the integral is evaluated over 
the whole field fl. Assuming a standard cloud composi- 
tion (63% hydrogen, 36% hehum, and 1% dust), we find 
H = 1.37 and a total mass M = (11 000 ± 2600) Mg. The 
error is mainly due to the uncertainty on the distance of 
the cloud. Note that in this analysis we only included mea- 
surements above b > -1-3° (see discussion above in Sect. 

121). 

Our mass estima te apparently compa res well with the 
independent one of lOnishi et"ari l)l999j) . M ~ 10^ Mq. 
No te, however, that if we use the cloud distance assumed 
bv lOnishi et al.l l|l999^ . d = 160 pc, we obtain a larger 
mass, M = 1.7 x 10* Mg, i.e., the CO derived mass is 
only about 65% the dust derived mass. As discussed in 
Sect. |S1 this discrepancy can in principle be attributed to 
(1) the insensitivity of ^^CO to low column densities, (2) 
to the saturation of ^^CO in the dense cores of the cloud, 
and (3 ) to a relatively small X-factor used bv lOnishi et alJ 
lll999l). We can rule o ut the latter as the X-factor used 
bv lOnishi etin (l1999^ (2.8 x 10^° cm'^ K"^ km"^ s) is 
virtually coincident (96%) with the one derived in this 
paper (2.91 x 10^° cm'^ K"^ km"^ s). 

To investigate the source of discrepancy in the mass 
estimates we present in Fig. |23the relationship between 
the integrated mass distribution and the extinction in Ak- 
We also plot, as a dashed fine, the smoothed and clipped 
extinction map that was compared to the CO data (see 
Fig. I21|l . Note that regions with extinction larger than 
Ak > 0.6 magnitudes (where the CO-dust correlation 
breaks at higher column densities) account for about 20% 
of the total mass (dashed line). Similarly, note how regions 
with Ak < 0.25 mag, the column density threshold below 
which CO is not sensitive to H2, account for about 15% 
of the cloud mass. The total fraction of dust mass that is 
missed by the CO is then about 35%, which, within the 
approximations, solves the discrepancy found between the 
dust and CO derived mass. 

7. Conclusions 

The main results of this paper can be summarized as fol- 
lows: 

— We used 4.5 million stars from the 2MASS point source 
catalog to construct a 8° x 6° extinction map of the 
poorly studied Pipe nebula, a molecular cloud complex 
seen in projection towards the galactic bulge. The map 
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Fig. 27. The cumulative mass enclosed in isoextinction 
contours. The plot has been constructed using only extinc- 
tion measurements at galactic latitude & > 3° (cf. Fig. 120(1 . 
The solid line shows the cumulative mass from the whole 
Nicer map, the dashed black line the cumulative mass 
from the downgraded (and clipped) map of Fig. 1211 right, 
and the dashed grey line the cumulative mass derived from 
the NANTEN observations (see Fig.EHright). For this lat- 
ter mass estimate we used the conversion formula (|15|l . 



has a resolution of 1 arcmin and has a Scr detection 
level of 0.5 visual magnitudes. 

— We combined the 2MASS data with Tycho and 
Hipparcos parallaxes to obtain a distance of 130±15 pc 
to this molecular complex. We estimated, for this dis- 
tance, a total mass of ^ lO'* Mq for the Pipe complex. 
This makes the Pipe nebula one the closest star form- 
ing region to the Earth of its type, closer than the 
Ophiuchus and the Taurus complexes. 

— We studied in details the statistical properties of the 
extinction map obtained through the Nicer method. 
We also con firmed the relation originally observed by 
iLada et all (^994) for IC 5146, suggesting that there is 
unresolved structure at the highest extinctions in our 
map. 

— We compared our near-infrared st udy with the CO 
observations of lOnishi et alJ l|l999l) . and derived fit- 
ting formulae to relate the ^^CO column density with 
the Ak extinction. We also derived the dust-to-^^CO 
ratio, and showed that if a normal infrared redden- 
ing law is assumed, then the derived X- factor is as 
large as 2.91 x 10^° cm^^ K^^ km~^ s in the range 
Ak e [0.1,0.6] mag. 

— We found that ^^CO is only sensitive to about 65% 
of the total dust mass. About half of the missing 
mass is not traced by CO at column densities be- 
low Ak < 0.25 mag, and half is not traced at column 
densities above Ak > 0.6 mag, where the line begins 
to saturate. There is an apparently irreducible uncer- 
tainty in the dust-CO correlation of about 1 magni- 
tude of visual extinction. This uncertainty seems to be 



independent of the cloud or the interstellar radiation 
field. 

— We took advantage of the large number of background 
sources to accurately measure the NIR reddening law 
for this cloud, and obtained E{J — H) — (1.85 ± 
Q.lb) E(H — K) , in very good agreement with the stan- 
dard lRieke fc L ebofskv (1985|) reddening law. 

— Finally from analysis of the JHK color-color diagram 
for the Pipe region we identify a large population of 
red stars whose colors are distinct from those of typical 
reddened background giants. These stars are spatially 
distributed over the entire observed field, preferentially 
located at lower Galactic latitudes and are not asso- 
ciated with the molecular cloud. These stars are the 
brightest stars detected in the field and have a very 
narrow distribution in magnitude space. These prop- 
erties are similar to those of Galactic OH/IR stars. 
Our observations thus appear to have provided one of 
the largest samples of such stars yet discovered. 
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